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Abstract. The spin-boson model, often used in NMR and ESR physics, quantum optics and spin- 
tronics, is considered in a solvable limit to model a spin one-half particle interacting with a bosonic 
thermal bath. By applying external pulses to a non-equilibrium initial state of the spin, work can be 
extracted from the thermalized bath. It occurs on the timescale 3? 2 inherent to transversal ('quan- 
tum') fluctuations. The work (partly) arises from heat given off by the surrounding bath, while the 
spin entropy remains constant during a pulse. This presents a violation of the Clausius inequality 
and the Thomson formulation of the second law (cycles cost work) for the two-level system. 



Introduction. After E.L. Hahn discovered the spin-echo in NMR physics [jp, it was 
soon considered to be a counterexample for irreversibility and for the second law [jT], 
0; for a recent discussion see e.g. [0, |3|]. In the present contribution we will show 
in a different context that NMR-physics contains quantum effects which should be 
interpreted as limits of the second law @|. Within the same general program we recently 
analyzed the thermodynamics of the Caldeira-Leggett model for a quantum harmonic 
oscillator coupled to a thermal bath |Q]. At low temperatures various formulations of 
the second law are violated: the Clausius inequality cti? < TdS is broken, the rates of 
energy dispersion and entropy production can be negative, and certain cycles are possible 
where heat extracted from the bath is fully converted into work ("perpetuum mobile"). 
The present analysis of the spin-boson model reveals a different quantum mechanism 
limiting the validity of the second law. 

The Hamiltonian of the problem reads: 

Jif — J4?(A) = J%? s + J4? B + (1) 
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This is a spin-^ interacting with a bath of harmonic oscillators (spin-boson model [|7|, |8p); 
J^, and stand for the Hamiltonians of the spin, the bath and their interaction, 
respectively. a x , o y and 6v = —id x 6 y are Pauli matrices, and and d k are the creation 
and annihilation operators of the bath oscillator with the index k, while the g k are the 
coupling constants. For an electron in a magnetic field B, £ = g}l B B is the energy, with g 
the gyro-magnetic factor and \i B the Bohr magneton. We shall consider a situation where 
A(f) = for almost all times. This is a prototype of a variety of physical systems [|7|], and 



known to be exactly solvable [|7], |8p, since the z-component of the spin is conserved, and 
with it the spin energy. Physically it means that we restrict ourselves to times much less 
than 3^ (relaxation time of 6 Z ). In ESR physics [|9|] the model represents an electron-spin 
interacting with a bath of phonons, for NMR it can represent a nuclear spin interacting 
with a spin bath, since in certain natural limits the latter can be mapped to the oscillator 
bath. In quantum optics it is suitable for describing a two-level atom interacting with a 
photonic bath JlOt]. 

Starting from general physical arguments [Q], one typically takes the quasi-Ohmic 
spectral density of the bath J(co) = T.kSk^i^k ~ a) )/ fi^k) = 8^ exp( — co/r) / %, where 
g is the dimensionless damping constant and the exponential cuts off the coupling at 
0) 3> r, the maximal frequency of the bath. As usual, the thermodynamic limit for the 
bath has been taken here. 

Since A = 0, one has conservation of o z (t) = (J z (0) (in the Heisenberg picture). Due 
to this one has 0]: 

^g k [al(t)+a k (t)] = f,(t)-a z G(t), (2) 
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where r) (t) is the quantum noise operator, and where the structure of G(t) shows that 
1 /r is the relaxation time of the bath. 

Separated initial state. To describe situations, where the spin was suddenly brought 
into the contact with the bath, e.g. an electron injected into semiconductor, atom injected 
into a cavity, or exciton created by external radiation, we make the assumption that 
initially, at t = 0, the spin and the bath are in a separated state, the latter being Gibbsian 
at temperature T = 1//3: p(0) = p s (0) ® exp(— f3j^ B )/Z B , where p 5 (0) is the initial 
density matrix of the spin. In this situation the quantum noise is stationary and Gaussian 
with average zero. The Heisenberg equation h~6 ± = i [±e + i)(t) —G(t)] 6± for <7 ± = 
& x ± / Cy can be solved exactly with the result [Q] : 



where (f)(t)fj(0) +f)(0)f)(?)) = 2ft 2 co Q = e/ti, and (o x (t)), (<7 y (t)) are determined 
via the real and imaginary part of (&±(t)). For t ^> l/T @ brings £(t) ~ t/ ^ = 
H/ (gT). £?2 can mus be identified with the transversal decay time. 

The density matrix of the spin reads P s = j + %Y,k=x yz^&W) ^ ts von Neumann 
entropy equals 5 vN = — trp 5 lnp 5 = —p } \np x —p 2 \np 2 , where p i2 = \ ±\\{a)\, d = 
{o x ,c y ,o z }. In the course of time \ (o(t)) \ decays to |(d z (0))|, which makes the von 
Neumann entropy increase. Since there is no heat flow - the energy is conserved - this 
is in agreement with a formulation of the second law: the entropy of closed system, or 
of an open system without energy transfer (the spin in contact with the bath), cannot 
decrease. 

A sudden pulse. So far we considered the Hamiltonian (HJ) with A = 0. A fast rotation 
around the x-axis is described by taking A ^ during a short time 8 l ; this is called a 



fast pulse [^]. If A ~ is large, the evolution operator describing the pulse becomes 
U l = exp(— id 1 J^ 3 (A) /H) ps exp(^' 0d x ) + ^(S^, where 9=—8 l A/H is the rotation angle, 
t/r dy Z t/j = dv sin ± 6y cos 0. During the sudden switchings of A(f ) from to A and 
from A to 0, the state of the system does not change, so p(t + 5j) = U 1 p(t) . The 
work done by the source is the change of the total energy, 

W 1 (t) = tr [p (t) (Jf(£) -J^)+p(t + 5 1 )(je-J^(A))]=\xp(t) {U ] ; l J^U l - Jf?) , 

since [U V J?(A)] = 0. 

Our main interest is work extraction from the bath. In order to ensure that the pulse 
does not change the energy of the spin, we first consider the case e = 0, where the spin 
has no energy. For small g, = — n/2 and t ^> 1/T the work appears to be 

= £ + (4) 

If for a fixed t, temperature is neither too large nor too small, 7e _f / ^ > HT/k, work 
can be extracted (W l < 0), provided the spin started in a coherent state (<7 X (0)) = — 1. 
This possibility to extract work from the equilibrated bath (tT ^> 1) contradicts to 
the Thomson's formulation of the second law out of equilibrium. It disappears on the 
timescale because then the spin looses its coherence, (o x , y (t)) — > 0. Notice that 
any combination of ±n pulses (this is a classical variation, since the coherence is not 
involved) can extract work only from a non-thermalized bath, i.e. for times ~ 1 /T. Thus, 
the effect is indeed essentially quantum. 

Initial preparation via a rotation. Our approach also allows to consider a specific, 
well controllable non-equilibrium initial state: a Gibbsian of the total system, p G — 
exp(— )/Z, in which at t = the spin is rotated ("zeroeth pulse") over an angle — \n 
around the y-axis, p(0) = U Q p G U l , with U Q = exp(—ind y /4). This maps 6 X — > c z , 
6 Z — > —d x . Such a state models the optical excitation of the spin, as it is done in NMR 
and spintronics. Though p(0) does not have the product form mentioned above (7), the 
problem remains exactly solvable. Taking = — \n one now gets 
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— sin co Q t + — cos co Q t 



tanh^e-^ (5) 



where y{t) = {g/n) arctanjT/] arises from friction, with 7(0°) = g/2. Typically g is 
small, so work is extracted (W 1 < 0) when the sine function is positive. The work 
decomposes as W 1 = AU — AQ, into the change in spin energy due to the pulse, AU ~ 
— (e/2) sina) f tanh[/3e/2] e^l'^-, and the heat absorbed from the bath 
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Notice its similarity with —W l of Eq. (|]). An interesting case is where work is performed 
by the total system (W l < 0) solely due to heat taken from the bath (AQ > 0, AU = 0). 
This process, possible by choosing t w 2nn/(O with integer n, can be considered as a 



cycle of a perpetuum mobile, forbidden by folklore minded formulations of the second 
law. Indeed, under a rotation the length \{o)\, and with it the von Neumann entropy, 
is left invariant, so one has a process with AQ > , A5 vN = 0, which violates the 
Clausius inequality AQ < TAS wN . The work needed at time zero to rotate the spin is 
W = (e/2) tanh[/3e/2] +gHT/ (2k), representing the work done on the spin and on the 
bath, respectively. It can be verified that the total work W + W l is always positive, so 
Thomson's formulation for a cyclic change [§| (here: the combination of the pulses at 
time t = and t) starting from equilibrium is obeyed. 

Two pulses in a rotated initial Gibbsian state. If there are many (very weakly inter- 
acting) spins, each in a slightly different external field, there appears an inhomogeneous 
broadening of the co = e/H line, for which we assume the distribution 



having average cb Q and inverse width 3~ 2 , typically much smaller than 3^. In this case 
the gain for a single pulse is washed out, leaving only the loss AQ = —gh~T /2k, so two 
pulses are needed. We consider again the rotated initial Gibbsian state, and perform a 
first — Ik pulse around the x-axis at time t l and a second jK pulse at time t 2 = t l + x 
(the standard \k,k combination would not expose an interesting role of the bath). In the 
regime of small g and large fj ~ 9" 2 the work in the second pulse is 

W 2 = ^-^i/^esina) Ttanh^ (8) 
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At moderate times only slowly oscillating terms survive. They are the ones that involve 
At = t 2 — 2t l . For the total work W = W 1 +W 2 the averaging over co Q brings 

K 4 2 

{ w n sin co n At + — - b \ ' ( 1 H n u ) cos co n At } 
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For At near 2Kn/(b Q such that the odd terms cancel, this again exhibits work extracted 
solely from the bath. 

Feasibility. Let us present several reasons favoring the feasibility of the proposed 
setups: 1) Two-level systems are widespread, because many quantum system act as 
two-level systems under proper conditions; 2) Detection in these systems is relatively 
easy, since already one-time quantities (o(t)) completely determine the state; 3) The 



harmonic oscillator bath is universal [TT|]; 4) Work and heat were measured in NMR 
experiments more than 35 years ago []12]; 5) Our main effects do survive the averaging 
over disordered ensembles of spins, thus allowing many-spin measurements. 6) The 
ongoing activity for implementation of quantum computers provides experimentally 
realized examples of two-level systems, which have sufficiently long 3T 2 times, and 



admit external variations on times smaller than 5^ : (/) for atoms in optical traps S? 2 ~ 1 s, 
1/r ~ 10~ 8 s, and there are efficient methods for creating non-equilibrium initial states 
and manipulating atoms by external laser pulses [jT3]; (ii) for an electronic spin injected 
or optically excited in a semiconductor ^ ~ 1 /is [14|]; (Hi) for an exciton created in a 
quantum dot ^ ~ 10~ 9 s [ JT5| ] (in cases (ii) and (Hi) 1/r ~ 10~ 13 s and femtosecond 
(10~ 15 s) laser pulses are available); (iv) in NMR physics ^ 2 ~ 10~ 6 — Is and the 
duration of pulses can be comparable with 1/r ~ 1 /is. 

In conclusion, we have analyzed for the spin-boson model the validity of some 
non-equilibrium formulations of the second law. The model is relevant for almost any 
branch of condensed matter, where two-level systems are described: NMR and ESR 
[§], Josephson junctions [0], quantum optics [JT0[] . Our main finding is that quantum 
coherence puts definite limits on the validity of non-equilibrium Thomson's formulation 
of the second law and on the validity of Clausius inequality. The effects disappear in the 
classical limit. The detailed discussion on the feasibility of the obtained effects can be 
found in [01 . 
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